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, ABSTRACT 

In the current ACDM cosmological scenario, TV-body simulations provide us with a 
Universal mass profile, and consequently a Universal equilibrium circular velocity of 
(*C) • the virialized objects, as galaxies. In this paper we obtain, by combining kinematical 

data of their inner regions with global observational properties, the Universal Rotation 
Curve (URC) of disk galaxies and the corresponding mass distribution out to their 
virial radius. This curve extends the results of Paper I, concerning the inner luminous 
i-G ■ regions of Sb-Im spirals, out to the edge of the galaxy halos. 
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£j I 1 INTRODUCTION 

(Z3 ■ 

. Rotation curves (hereafter RCs) of disk galaxies do not show any Keplerian fall-off and do not match the distribution of the 
^ ' stellar (plus gaseous) matter. As a most natural explanation, this implies an additional invisible mass component (Rubin et 
•i— ( ■ al. 1980; Bosma 1981, Persic & Salucci, 1988) that becomes progressively more conspicuous for the less luminous galaxies 
' (e.g.: Persic & Salucci 1988, 1990; Broeils 1992a). Moreover, the kinematical properties of Sb-Im spirals lead to the concept of 
^ . the Universal Rotation Curve (URC) implicit in Rubin 1985, pioneered in Persic and Salucci, 1991 and set in Persic, Salucci 
& Stel (1996, hereafter PSS, Paper I): RCs can be generally represented out to Ri, the outermost radius where data are 
available, by Vurc(R; P), i.e. by a universal function of radius, tuned by some galaxy property P. P can be a global property 
such as the luminosity and the disk or halo mass or a well defined local quantity like V op t- In any case it serves as the galaxy 
identifier. In PSS individual RCs and a number of coadded RCs proved the URC paradigm being well fitted by an analytical 
Curve, VuRc(r/R op t, ifl a function which is the sum in quadrature of two terms: Vurcd and Vurch, each representing the 
disk or halo contribution to the circular velocity: 

Vurc = Vurcd + Vurch (1) 
The stellar component was described by a Freeman disk (Freeman 1970) of surface density Sjj(r) = ^—§3 e~ r ^ RD and 
contributing to the circular velocity V as: 

V3rcd(x) = ±^^(3.2x) 2 (I K - hK x ) (2a) 

where x = r/R opt and I n and K n are the modified Bessel functions computed at 1.6 x. The dark matter component (with 
Vurch{ t ) = ° Mj ^ < ' ) was described by means of a simple halo velocity profile : 

Vurch{x) = ^-Gp a 2 x 2 /(a 2 + x 2 ) (2b) 

The above implies a density profile with an inner flat velocity core of size ~ aR op t, a central density po, an outer r~ 2 decline. 
The sum of the contributions (2a) and (2b) well fit all the PSS data with po, a 2 specific functions of luminosity (see PSS). 
bet us remind that disk masses Md of spirals were found in the range 10 9 M < M D < 2 x 10 11 M Q . 



1 The reader is directed to PSS for the details of the procedure. 

2 We define the "disk size" R op t = 3.2 Rd 
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The URC for the purpose of this work matches well the individual RCs of late type spirals (see also Appendix for 
a discussion). It is useful to express the URC paradigm in the following way: at any chosen radius, the URC predicts the 
circular velocity of a (late type) spiral of known luminosity and disk scale-length, within an error that is one order of magnitude 
smaller than the variations it shows i) at different radii and ii) at any radius, with respect to objects of different luminosity. 

Let us remind that the Universal curve built in PSS holds out to Ri, uses the luminosity as the galaxy identifier and the 
disk scale-length as a unit of measure for the radial coordinate. We will label it as URCo to indicate it as the first step of a 
definitive function of the dark radial coordinate, able to reproduce the observed RCs of spirals. URCo provides fundamental 
knowledge on the mass distribution in spirals, while it suffers from three limitations: 1) it strictly holds only in a region 
extended less than 5% the DM halo size (see below) 2) the velocity profile of the halo component, valid out to Ri, cannot be 
extrapolated to radii of cosmological interest 3) it identifies objects by their luminosities, rather than by their virial masses. 
Let us point out that the URCo has been often and successfully used as an observational benchmark for theories, but this, 
only for R < Ri and after that a relation between the halo mass and the galaxy luminosity was assumed. 

On the other side, high-resolution cosmological N-body simulations have shown that, within the (A) Cold Dark Matter 
(CDM) scenario, dark halos achieve a specific equilibrium density profile characterized by a universal shape and, in turn, an 
universal halo circular velocity (Navarro, Frenk & White, 1997, NFW), Vnfw(R, M v i r ) in which the virial mass M„; r and 
virial radius R vir are the galaxy identifier and radial coordinate. 

, . M vir c 2 g(c) , . 

Mr) = ^LWT^' (3a) 

where x = r/R v i r is the radial coordinate, c is the concentration parameter, and g(c) — [ln(l + c) — c/(l + c)] _1 . The parameter 
c is found to be a weak function of the halo mass, given by c « 14 (M„ lr /10 n M ) ' (Bullock et al, 2001, Dutton, 2006, 
Gnedin 2006). This leads to 

U N 2 P wM = K 2 lr -^-^ , (36) 
9(c) x 

with V vir — V(R v i r ). It is interesting to note that in this scenario the present-day circular velocity, which also includes a 
baryonic component arranged in a disk, is predicted to be a Universal function of radius, tuned by few galaxy parameters 
(Mo, Mao & White 1998). However, it is well known that observations of spiral galaxies favor density concentrations lower 
than those predicted for CDM by Eq. (3a): DM halos detected around spirals do not show the NFW central cusp in favor of 
a core- like structure (Gentile et 2005; van den Bosch & Swaters 2001; Swaters et al. 2003; Weldrake et al. 2003; Simon et al. 
2005; Donato et al, 2004; Gentile et al. 2007). 

Therefore, the reconstruction of the mass distribution of DM halos from observations in parallel with that emerging from 
N-body simulations is required not only as a normal scientific routine, but also in view of a theory-vs-observations likely 
disagreement. 

As an alternative to the simulation method, we will support the URC paradigm by means of a set of proper observational 
data and we will derive an analytical form for this curve, valid from the galaxy center out to its virial radius and characterized 
by the halo mass as the galaxy identifier. In detail, we extend/improve the URCo in PSS a) by adopting a different halo profile, 
proper to describe the halo distribution out to the virial radius, b) by using a number of RCs substantially more extended 
than those in PSS and c) by exploiting the relationship between the disk mass M D , and the virial galaxy mass M vir , recently 
obtained by Shankar et al. (2006). This will allow to build an "observational" Universal Curve, Vurc{R\ M v i r ), extended out 
to R v ir and having the virial mass as the galaxy identifier. This curve is the observational counterpart of the universal ACDM 
NFW N-body generated profile. 

While pointing that the concept behind the Universal Rotation Curve may be valid also for galaxies of different Hubble 
Types (see Salucci and Persic, 1997), but a number of issues are still open and will be dealt elsewhere: 

i) Sa galaxies amount, by number, to less than 10% of the whole spiral population, and are important objects in view of 
the dual nature of their stellar distribution. They show RC profiles with a clear systematics with luminosity (Rubin et al. , 
1985), but, not unexpectedly, with some difference from those of the URCo (Noordmeer, 2007). 

ii) Dwarf spirals with V op t < 50km/s are not well studied and included in the URC yet, also because in these objects the 
RCs do not coincide with the circular velocity, being significant the complex asymmetric drift correction. 

iii) The kinematical properties of spirals of very high stellar disk mass are not presently investigated with a suitably large 
sample. 

iv) A possible additional URC physical parameter (e.g. the surface stellar density) to take care of the (small) variance of 
the RCs profiles that seems to be unaccounted by the luminosity. 

Finally, let us remind that, in a flat cosmology with matter density parameter Qm = 0.27 and Hubble constant Ho = 
71 km s _1 Mpc -1 , at the present time, the halo virial radius Rvir-, i.e. the size of the virialized cosmological perturbation of 
mass M v i r scales with the latter as: 

^- 25 Kio^) 1/3kpc - w 
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Figure 1. The core radius vs. virial mass relations for the SE sample and the present work (solid line), 
(see e.g. Eke et al. 1996). 



2 THE UNIVERSAL HALO VELOCITY PROFILE 

We assume that the DM mass distribution is described by the Burkert (1995) profile 

P( r ) = 7 ^P-o 57 . ( 5ffl ) 

(r + ro) {r 2 + rl) 

where ro is the core radius and po the effective core density, in principle two independent parameters. Correspondingly, the 



ln(l + £) -tan- 1 (^) + i ln(l + ^)j with M = 1.6 po rg, so 



total halo mass inside radius r is given by Mh{t) = 4 Mo 
that: 

Vi RCH{ r) = 6.4 G ^%(l + -1) - tan- (I.) + ^„[l + (I.) 2 ] } . (56) 

Inside Ri this profile is indistinguishable from the halo term (2b) in the URCo (Salucci & Burkert, 2000, Gentile et al 2004). 
At larger radii, the mass diverges only logarithmically with radius and converges to the NFW velocity profile, provided that 

To <C Rvir- 

We fit the set of individual and coadded RCs of PSS with Vurc{R', Md, Po,ro) and derive the model parameters Md, po, 
ro (see Salucci and Burkert, 2001): 

lo S - -23-515 - 0.964 ( j,^-) (6a) 



10 11 M 



and 



r in-24 -2/3 -(rn/27) 2 -3 r r .r\ 

po — 5 x 10 r ' e y 'gem . (6b) 

Eqs. (l)-(2a)-(5b)-(6) define the URC out to Ri, Vurc{R, Md,vo), from the "baryonic perspective". Let us notice that, as 
result of the RC mass modeling, ro, differently from Md and po, has quite large fitting uncertainties, viz. oro/Vo = 0.3 — 0.5. 
Following our empirical approach, we do not extrapolate the URCH determined inside Ri out to R V i r 2> Ri ■ In that this will 
be uncertain besides of unknown validity. This quantity will be derived in the next section. 



3 THE URC OUT TO THE VIRIAL RADIUS 

We overcome the two main limitations of the URCo, its problematic extrapolation between Ri and R V i r and the uncertainty 
in the estimate of the core radius, by determining the latter by means of a new outer observational quantity, the halo virial 
velocity V v i r = [GM v i r /R v i r (M v i r )] 1 ^ 2 , related to the virial mass through Eq. (4). In detail, we obtain V v i r from the disk 
mass, suitably measured from inner kinematics through its relationship to the virial mass found by Shankar et al. (2006) 
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that is well represented (i.e. within 5%, the actual relation we use is shown in Fig. 2 and given in the Code indicated in the 
Discussion) by: 

M -oo Ml n 10 u [*W(3 IP" Me)] 3 ' 1 

This relationship is a consequence of the existence of i) the universal stellar mass function, (Bauldry et al 2004, Bell et al 
2003) and of ii) the cosmological halo mass function as indicated by iV-body ACDM simulations (f| Let us notice that this 
relationship is obtained without assuming any halo density profile, so that it can be combined with the mass modelling of the 
inner kinematics. 

Let us first derive RD(M vir ), the disk scale-length as a function of the halo mass, by inserting Eq. (7) in the relationship: 
log g = 0.633 + 0.379 log ^ + 0.069 (log ^) ' (8) 

obtained in PSS. We note that no result of this work is affected by the observational uncertainties on the relationship Eq. (8). 

It is worth to compute the radial extrapolation needed to reach R v i r from Ri = 6Rd, a quantity that can also be 
associated with the baryonic collapse factor F = R V i T jRo\ we find F ~ 90 — 15 log-^Yfj—, i-e. about 25 times the disk size 
~3R D . 

Eq. (7) in combination with Eq. (4) allows us to add, to the PSS set of kinematical data leading to the URCo, a new 
observational quantity: V v i r (Mo) = GM v i r / R v i r , relative to the virial radius. Then, we determine the core radius not from 
the inner kinematics, but as the value of ro for which the velocity model described by Eqs. (1), (2a), (5b) and (6a), matches 
(at Rvir) the virial velocity V v i r given by Eq. (7) and (4). Let us write this as: 

□ — 7TT \ = VuRCH[Rvir(M vir ); p (M vir ) , r )] (9a) 

where p(M v i r ) as a short form for p(MD(Mmr)) with Eq. (7) inserted in Eq. (6a). From the PSS inner kinematics we get the 
values of po and Md according to Eq. (6a), so that Eq. (9a) becomes an implicit relation between ro and Mvir (ci,C2 are 
known numerical constants): 

M vir = Cl ^^ (Jl + - tan- + hull + f^Y] } . (96) 

oiM ' ^ \ r J \ r ) 2 L V r / J J 

The above can be numerically solved for any M v %r, and the solution can be approximated by: 

log (ro/kpc) ~ 0.66 + 0.58 log (M^ ir /10 X1 M ) . (10) 

(an higher order approximation is given in the Code indicated in the Discussion). Let us stress that the present derivation 
of ro is very solid with respect to observational uncertainties: errors up to a factor 2 in M v i r in eq (7) trigger errors in ro 
lower than 40%, and errors in the outer halo velocity slope (0.1 < R/R V i r < 1) lower than 0.1. This is certainly smaller than 
the scatter of values with which this quantity is found by iV-Body simulations and by SPH/semi analytical studies of galaxy 
formation including the baryonic components. 

It is worth investigating a number of recently published super-extended (SE) RCs (Donato et al. 2005, Gentile et al 2004, 
Salucci et al 2003). They reach a radius larger than 5% (and up to 15%) of the virial radius, i.e. a radius at least twice as 
extended as those of the sinthetic curves in PSS. The mass modelling of these SE RCs (made in the original papers) shows 
an ro vs. M v i r relationship that is in good agreement with Eq. (10). Relation (10) and the above individual values differ by 
20% -40% from those determined from the inner kinematics alone and given by means of eq (6b). Since in this paper (also 
because Ri ~ ro), eq (6) is considered a prediction of the inner mass modelling rather than an actual measurement, such good 
agreement indicates the soundness of the PSS mass modelling. 

Let us notice that only for a range of values of the crucial quantity V(Ri)- Vvir, with the first term obtained by the inner 
kinematics and the second one via Eqs. (4) and (7), there is a solution for eq.(9b), therefore, the existence of eq (10) and the 
agreement of the values of the Burkert core radii, measured independently at 0.05JiW, Q.lR v i r and R v i r are important tests 
passed by this profile. 

Then, by means of Eqs. (1), (2a)-(5b)-(6a)-(10), we construct the full URC, extended out to the virial radius and with 
the virial mass as the galaxy indicator. It is useful to show the relationships we use (see Fig. 2). The mass model includes a 
Burkert DM halo of central density po, of core radius of size ro and a Freeman disk of mass Md- The URC fits nicely the 
available velocity data out to Ri and it is valid out to the virial radius, where it exactly matches V v ir- Moreover, since M v i r 
is the quantity that in theoretical studies identifies a galaxy, we overcome the main limitation of URCo- 

We consider all of the three coordinate systems r,r / Ro,r / R v i r equivalent to represent the main structural properties of 



3 There is no inconsistency in adopting the ACDM halo mass function and cored halo mass models, in that the latter can be formed 
astrophysical from the cosmological cuspy ones. Since both functions account the same cosmological objects, the Jacobian of their 
transformation defines a relation between the disk and virial mass in spirals (see Shankar et al. 2006). 
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Figure 2. The various relationships used in this paper. 
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Figure 3. The Universal Rotation Curve in physical units. Each curve corresponds to M„; r = 10 11 lO n / 5 M0, with n = 1 . . .9 from the 
lowest to the highest curve. 

the mass distribution in spirals, but each of them showing some particular aspects. More specifically, it is then possible and 
useful to build several "URCs" i.e. VuRc(r/R CO o', P), where P is a galaxy identifier (M D , M vir , L) and R coo a radial coordinate 
(r,r/R D ,r/R v i r ). Although not all these URCs are independent in a statistical sense, they all are relevant in that they all 
well reproduce the individual RCs and each of them highlights particular properties of the mass distribution. 

In Fig. 3 we show Vunc(r; M vir ), the URC in physical units with the objects identified by the halo virial mass; each line 
refers to a given halo mass in the range 10 11 M© <M v i r < 10 13 M ; the halo mass determines both the amplitude and the 
shape of the curve. Note the contribution of the baryonic component, negligible for small masses but increasingly important 
in the larger structures, mirrors the behavior of the M V i r — Md relation. The general existence of an inner peak is evident but, 
especially at low masses, it is due to both dark and stellar components. Remarkably, the maximum value of the circular velocity 
occurs at about 15 ± 3 kpc, independent of the galaxy mass: this seems to be a main kinematic imprint of the DM - luminous 
mass interaction occurring in spirals. Furthermore, Fig. 3 shows that the "Cosmic Conspiracy" paradigm has no observational 
support: there is no fine tuning between the dark and the stellar structural parameters to produce the same particular RC 
profile in all objects (e.g. a fiat one). Conversely, a number of relationships between the various structural parameters produce 
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Figure 4. The Universal Rotation Curve, normalized at its virial value VjjRC {Rvir), as a function of normalized dark radius x = R/R v { r . 
Each curve, from the highest to the lowest, corresponds to M v i r defined as in Fig. 3. The bold line is the NFW velocity profile (see text). 

a variety of RC profiles. Moreover, the peak velocity of the stellar component Vj^* = Vd(2.2.Rd) = GMd/Rd k, with 
k = const, is not a constant fraction of the virial velocity as is found in ellipticals, (i.e a oc Vvir), but it ranges between the 
values 1 and 2 depending on the halo mass. 

Moreover, as in the NFW (and Burkert) RC profiles, the URC profiles are found (moderately) decreasing over most of 
the halo radial extent. The paradigm of flat rotation curves is obviously incorrect even/especially intended as an asymptotic 
behavior at large radii. In fact, we find that both U(0.05iiW),the velocity at the farthest radius with available kinematical in 
PSS and V(3Rd), a main reference velocity of the luminous regions of spirals, are significantly (10% — 30%) higher than the 
(observational) value of V V % T - This rules out a V = constant extrapolation of the inner RCs out to regions non mapped by 
the kinematics and DM dominated regions. We note that this result is independent of the adopted halo density profile and is 
far from being granted on theoretical grounds. 

In Fig. 4 we frame the URC from a full DM perspective by plotting Vjjrc{R/ Rvir', M„j r ). We set the virial mass M v i r 
as the galaxy identifier and R/R V i r as the radial "dark" coordinate, thus normalizing the amplitudes by V V i r oc M^. This 
ensemble of curves, a main goal of the present work, is parallel to those emerging in iV-body simulations and aims to represent 
the actual velocity profiles of spirals. In these variables the DM halos are self-similar; the whole system is self-similar in 
the outer regions, while in the innermost 30% of the halo size the baryons have influenced the dynamics and broken the 
self-similarity. In these coordinates it easily emerges that the maximum of the RC occurs at very different radii, viz. at ~ 2Rd 
for the most massive objects and at ~ 10-Rd for the least massive ones. Then, no reference circular velocities, to be considered 
as the actual physical counterparts of the empirical velocities of the Tully-Fisher relationship, exist in actual galaxy rotation 
curves. 

In Fig. 5 we zoom into the URC to look for the inner (luminous) regions of spirals from a baryonic perspective: the URC 
is so expressed as a function of the "baryonic" radial coordinate r/Ro- This figure corresponds to Fig. 4 of PSS, with the 
important difference that here the virial mass, rather than the galaxy luminosity, is the galaxy identifier. Plainly, an inverse 
correlation between the average steepness of the RC slope and the halo mass holds, similar to the slope-luminosity relationship 
found by Persic & Salucci (1988). In this coordinate the stellar matter is closely self-similar, and the different shapes of the 
RCs curves are mainly due to the M v i r — Md relation. In the space defined by normalized circular velocity - dark radius - 
halo mass, spirals do not occupy random positions, but a well defined plane of very small thickness. We clearly see that, by 
filling only less than 10 _3 % of the available volume, the available kinematics of spiral galaxies defines the Universal Rotation 
Curve. Let us notice that, in principle, theories of the formation of spirals do not trivially imply the existence of such a surface 
that underlies the occurrence of a strong dark - luminous coupling. 

We now show the URC DM density distribution. In Fig. 6 we show it as a function of x and M V i r . For x < 10 _1 
the well known core-cusp discrepancy emerges, i.e. the DM density of actual halos around spirals is about one order of 
magnitude smaller and radially much more constant than the NFW predictions. At x > 0.4, for a concentration parameter 
c = 13(M„ir/10 12 Mq)~ ' 13 , the observed halo densities are consistent with the NFW predictions for halos of the same virial 
mass. Note that this is a direct test: for halos with density profiles at x > 0.5 very different from the Burkert or the NFW 
profiles, Eq. 9b does not have solution. 

More specifically, let us constrain the analytical form of the outer DM distribution. For 2 < r/ro < 18, the following 
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Figure 6. The URC halo density vs. the NFW halo density of objects of the same mass, as a function of normalized radius and virial 
mass. The axes labels are x, logM v i r /Mq and log(p/(g cm -3 )) 
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approximation for the Burkert and NFW profile holds (y = r/ro, e = 0) 

V URC «{y) = Wh(3.24) i 59+ ^ 119+6 . (11) 

Let us suppose that the actual outer DM velocity profile is different from the Burkert/NFW given by Eq. (11), i.e. t / 0. 
Then in Fig 7 we show that, even assuming large uncertainties in V V i r , in order to match both V(Ri) and Vvir, we must have 
e < 0.1. This is a first direct support for the Burkert and the NFW density law to be able to represent the outer regions 
(0.3 < r/R vir < 1) of DM galaxy halos. 

Notice that weak-lensing shear fields, at several hundreds kpc from the galaxy centers, are found compatible with the 
predictions of the NFW density profile, but cannot exclude non-NFW profiles ( Kleinheinrich ., et al., 2006 and references 
therein). 



4 DISCUSSION AND CONCLUSIONS 



In this paper, we have built the Universal Rotation Curve of spiral galaxies by means of kinematical and photometric data. 
We physically extended the URC, established for the inner region of galaxies in PSS out to R V i r and have been able to employ 
the virial mass M vir as the parameter that characterizes spiral galaxies, and the virial radius R V i r as a unit of measure for 
the radial coordinate. This URC is meant to be the observational counterpart of the NFW rotation curve, emerging from 
cosmological simulations performed in the CDM scenario. The URC yields the gravitational potential at any radius and it 
allows to link the local properties in the inner luminous regions with the global properties of the DM halos. 

DM halos have one (and likely just one) characteristic length scale, ro oc M®i®, which it is not naturally present in current 
scenarios of galaxy formation. Thus, they do not show any sign of an inner cuspy region of size r s oc M®it- The halo velocity 
contribution Vurch rises with radius like a solid body at r ~ 0, decelerates to reach a maximum at 3.24 ro from where it 
start to slowly decrease out to R V i r with a slope that it is consistent with that of the NFW halos. The main significance of 
the URC concerns the full mass distribution (MD). First, it is possible to immediately exclude the following scenarios (and 
combinations of them): i) individual behavior, every object has its own MD; ii) unique behavior, every object has almost the 
same MD. Instead, the MD in spirals shows a remarkable mass-dependent systematics: both the dark and the stellar matter 
are distributed according to profiles that are functions of the total mass M v i r (see Fig. 8). Finally, the DM halo becomes the 
dominant mass component in galaxies at different radii, according to the galaxy mass: from ~ lQ~ 2 R v i r for the lowest masses, 
to ~ 10~ l Rvi r for the highest ones. 

We write the equilibrium velocity of the halos around spirals as the following approximation of the relations in the 
previous Sect.: 



Vurci 



"Win 



1 + j(M vir )x 



tan 



j(M v 



1, 

+ -ln 



1 + j(M vir )x' 



with A(M viT ) = 0.406+1.08 /o 3 [M„i r /(lO n M )]^O.688 {log[M vir /(lO n M )]} 2 )+O.766 {log[M vir /(lO n M )]} 3 a,ndy(M vir ) = 



26.78 [M„ ir /(10 n M o )]- 



This is the observational counterparts of JV-body outcomes. 



A Mathematica code for the figures in this paper is available at: http://www.novicosmo.org/salucci.asp 



5 APPENDIX 

In this Appendix we discuss the observational evidence for the URC claim, the nature and the implications of which it is 
worth to clarify. The paradigm states that, when binned by luminosity, the RCs form a set of smooth, low-scatter synthetic 
curves, whose profiles and amplitudes are strong functions of the luminosity bin.0. Furthermore, the URC paradigm implies: 
i) rotation velocity slopes vs. rotation velocity amplitudes relationships (see fig (2) and (3) of PSS) and ii) a set of relations 
(Radial Tully Fisher relationship) holding at different radii x, defined as: 

log V(x) = a x M + b x 

where x = R/Rd and a x and b x are the fitting parameters, and M is the galaxy magnitude (Persic and Salucci 1991). 

Evidence for the URC claim and/or its above implications comes from: a) detailed analyzes of independent samples: 
Catinella et al 2006, (2200 RCs, see their Fig. 12), Swaters 1999, (60 extended RCs, see Chapter 4); b) independent analyzes 
of the PSS sample: Rhee (1996), Roscoe (1999); c) the finding of a very tight RTF in PSS and other three different samples 
Willick (1997, see below), Yegorova et al. (2007). 

The claim has been also tested by comparing the RCs of two samples of spirals (Courteau, 1998, 131 objects; Verheijen, 



4 The analytical form of the URC is built by assuming reasonable disk-halo velocity profiles, with three free parameters (V(Ropt), a, f3) 
that are obtained by x 2 -htting the sinthetic curves. 
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Figure 7. Halo velocities at Ri and R v i r (filled circles) vs. the URC-halo velocity, given by Eqs. (11) (solid line) and vs. velocity profiles 
with average logarithmic slope steeper or shallower by an amount e = 0.1 (dashed line). 
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Figure 8. The dark halo (solid line), the disk (short dashed line), and the total (long dashed) URC mass profile for reference masses of 
1O 11 M0 and 1O 13 M0. The distribution for the intermediate halo masses can be derived from section 2 

1997, 30 objects) with the circular velocities predicted by the URCo, once that the values of galaxy luminosity and disk length- 
scale are inserted in it. The face- value result of the test: 2/3 of the RCs are in pretty good agreement with the universal 
curve, while 1/3 show some disagreement, indicates that the URC is a useful tool to investigate the systematics of the mass 
distribution in spirals, but also it questions about its universality. However, while some of this disagreement may reflect an 
inefficiency of the URCo to reproduce the RCs, the actual performance of the URC is better than it is claimed. In fact, 
spurious data vs predictions disagreements are created in performing this test and precisely when they insert in URCo the 
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values of Lb and Rd, affected by (occasionally large) observational errors. By taking into account this effect the URCo success 
rate reaches 80% and more. 

Willick (f998) found, by studying a large sample of RCs, a radial variation of the scatter of the inverse RTF defined 
above and he interpreted it as an evidence against the URC. Let us show that this argument is incorrect and that au contraire 
the properties of the RTF support the URC paradigm. The increase/decrease of the scatter found is very small (Willick, 
1998): the scatter ranges from 0.065 dex (at 2Rd ) to 0.080 dex (at 0.5Rd and at 3Rd) and it implies, if totally intrinsic, a 
prediction error in log V{x) of (0.08 2 - 0.065 2 ) 5 = 0.04 dex. Moreover, some of the scatter increase/decrease is due to the 
larger random observational errors present in the outermost measurements; in fact, a refined analysis of the issue (Yegorova 
et al 2007) finds a smaller predicting error for three large sample of spirals. Therefore, from the RTF we have that, in the 
region considered, the luminosity statistically predicts the circular velocity at any radius and in any galaxy within an error of 
5% - 10% , a quantity much smaller than the variations of the latter in each galaxy and among galaxies. 
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, ABSTRACT 

In the current ACDM cosmological scenario, TV-body simulations provide us with a 
Universal mass profile, and consequently a Universal equilibrium circular velocity of 
(*C) • the virialized objects, as galaxies. In this paper we obtain, by combining kinematical 

data of their inner regions with global observational properties, the Universal Rotation 
Curve (URC) of disk galaxies and the corresponding mass distribution out to their 
virial radius. This curve extends the results of Paper I, concerning the inner luminous 
i-G ■ regions of Sb-Im spirals, out to the edge of the galaxy halos. 

Oh 1 

6 

£j I 1 INTRODUCTION 

(Z3 ■ 

. Rotation curves (hereafter RCs) of disk galaxies do not show any Keplerian fall-off and do not match the distribution of the 
^ ' stellar (plus gaseous) matter. As a most natural explanation, this implies an additional invisible mass component (Rubin et 
•i— ( ■ al. 1980; Bosma 1981, Persic & Salucci, 1988) that becomes progressively more conspicuous for the less luminous galaxies 
' (e.g.: Persic & Salucci 1988, 1990; Broeils 1992a). Moreover, the kinematical properties of Sb-Im spirals lead to the concept of 
^ . the Universal Rotation Curve (URC) implicit in Rubin 1985, pioneered in Persic and Salucci, 1991 and set in Persic, Salucci 
& Stel (1996, hereafter PSS, Paper I): RCs can be generally represented out to Ri, the outermost radius where data are 
available, by Vurc(R; P), i.e. by a universal function of radius, tuned by some galaxy property P. P can be a global property 
such as the luminosity and the disk or halo mass or a well defined local quantity like V op t- In any case it serves as the galaxy 
identifier. In PSS individual RCs and a number of coadded RCs proved the URC paradigm being well fitted by an analytical 
Curve, VuRc(r/R op t, L) , a function which is the sum in quadrature of two terms: Vurcd and Vurch, each representing the 
disk or halo contribution to the circular velocity: 

Vurc = Vurcd + Vurch (1) 
The stellar component was described by a Freeman disk (Freeman 1970) of surface density Sjj(r) = ^—§3 e~ r ^ RD and 
contributing to the circular velocity V as: 

V3rcd(x) = ±^^(3.2x) 2 (I K - hK x ) (2a) 

where x = r/R opt 2 and I n and K n are the modified Bessel functions computed at 1.6 x. The dark matter component (with 
Vurch(i~) = ° Mj ^ < ' ) was described by means of a simple halo velocity profile : 

Vurch{x) = ^-Gp a 2 x 2 /(a 2 + x 2 ) (2b) 

The above implies a density profile with an inner flat velocity core of size ~ aR op t, a central density po, an outer r~ 2 decline. 
The sum of the contributions (2a) and (2b) well fit all the PSS data with po, a 2 specific functions of luminosity (see PSS). 
bet us remind that disk masses Md of spirals were found in the range 10 9 M < M D < 2 x 10 11 M Q . 



1 The reader is directed to PSS for the details of the procedure. 

2 We define the "disk size" R op t = 3.2 R^ 
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The URC for the purpose of this work matches well the individual RCs of late type spirals (see also Appendix for 
a discussion). It is useful to express the URC paradigm in the following way: at any chosen radius, the URC predicts the 
circular velocity of a (late type) spiral of known luminosity and disk scale-length, within an error that is one order of magnitude 
smaller than the variations it shows i) at different radii and ii) at any radius, with respect to objects of different luminosity. 

Let us remind that the Universal curve built in PSS holds out to Ri, uses the luminosity as the galaxy identifier and the 
disk scale-length as a unit of measure for the radial coordinate. We will label it as URCo to indicate it as the first step of a 
definitive function of the dark radial coordinate, able to reproduce the observed RCs of spirals. URCo provides fundamental 
knowledge on the mass distribution in spirals, while it suffers from three limitations: 1) it strictly holds only in a region 
extended less than 5% the DM halo size (see below) 2) the velocity profile of the halo component, valid out to Ri, cannot be 
extrapolated to radii of cosmological interest 3) it identifies objects by their luminosities, rather than by their virial masses. 
Let us point out that the URCo has been often and successfully used as an observational benchmark for theories, but this, 
only for R < Ri and after that a relation between the halo mass and the galaxy luminosity was assumed. 

On the other side, high-resolution cosmological N-body simulations have shown that, within the (A) Cold Dark Matter 
(CDM) scenario, dark halos achieve a specific equilibrium density profile characterized by a universal shape and, in turn, an 
universal halo circular velocity (Navarro, Frenk & White, 1997, NFW), Vnfw(R, M v i r ) in which the virial mass M„; r and 
virial radius R vir are the galaxy identifier and radial coordinate. 

, . M vir c 2 g(c) , . 

Mr) = ^LWT^' (3a) 

where x = r/R v i r is the radial coordinate, c is the concentration parameter, and g(c) — [ln(l + c) — c/(l + c)] _1 . The parameter 
c is found to be a weak function of the halo mass, given by c « 14 (M„ lr /10 n M ) ' (Bullock et al, 2001, Dutton, 2006, 
Gnedin 2006). This leads to 

U N 2 P wM = K 2 lr -^-^ , (36) 
9(c) x 

with V vir — V(R v i r ). It is interesting to note that in this scenario the present-day circular velocity, which also includes a 
baryonic component arranged in a disk, is predicted to be a Universal function of radius, tuned by few galaxy parameters 
(Mo, Mao & White 1998). However, it is well known that observations of spiral galaxies favor density concentrations lower 
than those predicted for CDM by Eq. (3a): DM halos detected around spirals do not show the NFW central cusp in favor of 
a core- like structure (Gentile et 2005; van den Bosch & Swaters 2001; Swaters et al. 2003; Weldrake et al. 2003; Simon et al. 
2005; Donato et al, 2004; Gentile et al. 2007). 

Therefore, the reconstruction of the mass distribution of DM halos from observations in parallel with that emerging from 
N-body simulations is required not only as a normal scientific routine, but also in view of a theory-vs-observations likely 
disagreement. 

As an alternative to the simulation method, we will support the URC paradigm by means of a set of proper observational 
data and we will derive an analytical form for this curve, valid from the galaxy center out to its virial radius and characterized 
by the halo mass as the galaxy identifier. In detail, we extend/improve the URCo in PSS a) by adopting a different halo profile, 
proper to describe the halo distribution out to the virial radius, b) by using a number of RCs substantially more extended 
than those in PSS and c) by exploiting the relationship between the disk mass M D , and the virial galaxy mass M vir , recently 
obtained by Shankar et al. (2006). This will allow to build an "observational" Universal Curve, Vurc{R\ M v i r ), extended out 
to R v ir and having the virial mass as the galaxy identifier. This curve is the observational counterpart of the universal ACDM 
NFW N-body generated profile. 

While pointing that the concept behind the Universal Rotation Curve may be valid also for galaxies of different Hubble 
Types (see Salucci and Persic, 1997), but a number of issues are still open and will be dealt elsewhere: 

i) Sa galaxies amount, by number, to less than 10% of the whole spiral population, and are important objects in view of 
the dual nature of their stellar distribution. They show RC profiles with a clear systematics with luminosity (Rubin et al. , 
1985), but, not unexpectedly, with some difference from those of the URCo (Noordmeer, 2007). 

ii) Dwarf spirals with V op t < 50km/s are not well studied and included in the URC yet, also because in these objects the 
RCs do not coincide with the circular velocity, being significant the complex asymmetric drift correction. 

iii) The kinematical properties of spirals of very high stellar disk mass are not presently investigated with a suitably large 
sample. 

iv) A possible additional URC physical parameter (e.g. the surface stellar density) to take care of the (small) variance of 
the RCs profiles that seems to be unaccounted by the luminosity. 

Finally, let us remind that, in a flat cosmology with matter density parameter Qm = 0.27 and Hubble constant Ho = 
71 km s _1 Mpc -1 , at the present time, the halo virial radius Rvir-, i.e. the size of the virialized cosmological perturbation of 
mass M v i r scales with the latter as: 

^- 25 Kio^) 1/3kpc - w 
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Figure 1. The core radius vs. virial mass relations for the SE sample and the present work (solid line), 
(see e.g. Eke et al. 1996). 



2 THE UNIVERSAL HALO VELOCITY PROFILE 

We assume that the DM mass distribution is described by the Burkert (1995) profile 

P( r ) = 7 ^P-o 57 . ( 5ffl ) 

(r + ro) {r 2 + rl) 

where ro is the core radius and po the effective core density, in principle two independent parameters. Correspondingly, the 



ln(l + £) -tan- 1 (^) + i ln(l + ^)j with M = 1.6 po rg, so 



total halo mass inside radius r is given by Mh{t) = 4 Mo 
that: 

Vi RCH{ r) = 6.4 G ^%(l + -1) - tan- (I.) + ^„[l + (I.) 2 ] } . (56) 

Inside Ri this profile is indistinguishable from the halo term (2b) in the URCo (Salucci & Burkert, 2000, Gentile et al 2004). 
At larger radii, the mass diverges only logarithmically with radius and converges to the NFW velocity profile, provided that 

To <C Rvir- 

We fit the set of individual and coadded RCs of PSS with Vurc{R', Md, Po,ro) and derive the model parameters Md, po, 
ro (see Salucci and Burkert, 2001): 

lo S - -23-515 - 0.964 ( j,^-) (6a) 



10 11 M 



and 



r in-24 -2/3 -(rn/27) 2 -3 r r .r\ 

po — 5 x 10 r ' e y 'gem . (6b) 

Eqs. (l)-(2a)-(5b)-(6) define the URC out to Ri, Vurc{R, Md,vo), from the "baryonic perspective". Let us notice that, as 
result of the RC mass modeling, ro, differently from Md and po, has quite large fitting uncertainties, viz. oro/Vo = 0.3 — 0.5. 
Following our empirical approach, we do not extrapolate the URCH determined inside Ri out to R V i r 2> Ri ■ In that this will 
be uncertain besides of unknown validity. This quantity will be derived in the next section. 



3 THE URC OUT TO THE VIRIAL RADIUS 

We overcome the two main limitations of the URCo, its problematic extrapolation between Ri and R V i r and the uncertainty 
in the estimate of the core radius, by determining the latter by means of a new outer observational quantity, the halo virial 
velocity V v i r = [GM v i r /R v i r (M v i r )] 1 ^ 2 , related to the virial mass through Eq. (4). In detail, we obtain V v i r from the disk 
mass, suitably measured from inner kinematics through its relationship to the virial mass found by Shankar et al. (2006) 
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that is well represented (i.e. within 5%, the actual relation we use is shown in Fig. 2 and given in the Code indicated in the 
Discussion) by: 

Afo_ 2.3X10 M 01 + [AW(3 1Q11M0)] , 2 • (7) 

This relationship is a consequence of the existence of i) the universal stellar mass function, (Bauldry et al 2004, Bell et al 
2003) and of ii) the cosmological halo mass function as indicated by A^-body ACDM simulations . 3 Let us notice that this 
relationship is obtained without assuming any halo density profile, so that it can be combined with the mass modelling of the 
inner kinematics. 

Let us first derive R D (M vir ), the disk scale-length as a function of the halo mass, by inserting Eq. (7) in the relationship: 
log g = 0.633 + 0.379 log + 0.069 (log ' (8) 

obtained in PSS. We note that no result of this work is affected by the observational uncertainties on the relationship Eq. (8). 

It is worth to compute the radial extrapolation needed to reach R vir from Ri = 6Rd, a quantity that can also be 
associated with the baryonic collapse factor F — R vir /RD\ we find F m 90 — 15 W 10 iim q ' *' e ' arjout 25 times the disk size 
~ 3Rd- 

Eq. (7) in combination with Eq. (4) allows us to add, to the PSS set of kinematical data leading to the URCo, a new 
observational quantity: V v i r (Mo) = GM v i r / R v i r , relative to the virial radius. Then, we determine the core radius not from 
the inner kinematics, but as the value of ro for which the velocity model described by Eqs. (1), (2a), (5b) and (6a), matches 
(at R v i r ) the virial velocity V v i r given by Eq. (7) and (4). Let us write this as: 

GA1 

VuRCH[Rvir(M vir ); p (M vir ),r )] (9a) 



Rvir ( ^Ivi 

where p(M V ir) as a short form for p(MD{M V ir)) with Eq. (7) inserted in Eq. (6a). From the PSS inner kinematics we get the 
values of po and Md according to Eq. (6a), so that Eq. (9a) becomes an implicit relation between ro and M V i r (ci,C2 are 
known numerical constants): 

M v ir = Cl ^#(/n(l + ^) - tan- + hn\l + f^lY] } . (») 

coM ' I- V ro J V r ) 2 L V r / J J 

The above can be numerically solved for any M v i r , and the solution can be approximated by: 

log (ro/kpc) ~ 0.66 + 0.58 log (M^/10 11 M ) . (10) 

(an higher order approximation is given in the Code indicated in the Discussion). Let us stress that the present derivation 
of ro is very solid with respect to observational uncertainties: errors up to a factor 2 in M V i T in eq (7) trigger errors in ro 
lower than 40%, and errors in the outer halo velocity slope (0.1 < R/R V i r < 1) lower than 0.1. This is certainly smaller than 
the scatter of values with which this quantity is found by W-Body simulations and by SPH/semi analytical studies of galaxy 
formation including the baryonic components. 

It is worth investigating a number of recently published super-extended (SE) RCs (Donato et al. 2005, Gentile et al 2004, 
Salucci et al 2003). They reach a radius larger than 5% (and up to 15%) of the virial radius, i.e. a radius at least twice as 
extended as those of the sinthetic curves in PSS. The mass modelling of these SE RCs (made in the original papers) shows 
an ro vs. M v i r relationship that is in good agreement with Eq. (10). Relation (10) and the above individual values differ by 
20% -40% from those determined from the inner kinematics alone and given by means of eq (6b). Since in this paper (also 
because Ri ~ ro), eq (6) is considered a prediction of the inner mass modelling rather than an actual measurement, such good 
agreement indicates the soundness of the PSS mass modelling. 

Let us notice that only for a range of values of the crucial quantity V(Ri)- V V ir, with the first term obtained by the inner 
kinematics and the second one via Eqs. (4) and (7), there is a solution for eq.(9b), therefore, the existence of eq (10) and the 
agreement of the values of the Burkert core radii, measured independently at 0.05R v i r , Q.lRvir and R V i r are important tests 
passed by this profile. 

Then, by means of Eqs. (1), (2a)-(5b)-(6a)-(10), we construct the full URC, extended out to the virial radius and with 
the virial mass as the galaxy indicator. It is useful to show the relationships we use (see Fig. 2) . The mass model includes a 
Burkert DM halo of central density po, of core radius of size ro and a Freeman disk of mass Md- The URC fits nicely the 
available velocity data out to Ri and it is valid out to the virial radius, where it exactly matches V V ir- Moreover, since M v i r 
is the quantity that in theoretical studies identifies a galaxy, we overcome the main limitation of URCo- 

We consider all of the three coordinate systems r,r /Ro,r / R v i r equivalent to represent the main structural properties of 



3 There is no inconsistency in adopting the ACDM halo mass function and cored halo mass models, in that the latter can be formed 
astrophysical from the cosmological cuspy ones. Since both functions account the same cosmological objects, the Jacobian of their 
transformation defines a relation between the disk and virial mass in spirals (see Shankar et al. 2006). 
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Figure 2. The various relationships used in this paper. 
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Figure 3. The Universal Rotation Curve in physical units. Each curve corresponds to M„; r = 10 11 lO n / 5 M0, with n = 1 . . .9 from the 
lowest to the highest curve. 

the mass distribution in spirals, but each of them showing some particular aspects. More specifically, it is then possible and 
useful to build several "URCs" i.e. VuRc(r/R CO o', P), where P is a galaxy identifier (M D , M vir , L) and R coo a radial coordinate 
(r,r/R D ,r/R v i r ). Although not all these URCs are independent in a statistical sense, they all are relevant in that they all 
well reproduce the individual RCs and each of them highlights particular properties of the mass distribution. 

In Fig. 3 we show Vunc(r; M vir ), the URC in physical units with the objects identified by the halo virial mass; each line 
refers to a given halo mass in the range 10 11 M© <M v i r < 10 13 M ; the halo mass determines both the amplitude and the 
shape of the curve. Note the contribution of the baryonic component, negligible for small masses but increasingly important 
in the larger structures, mirrors the behavior of the M V i r — Md relation. The general existence of an inner peak is evident but, 
especially at low masses, it is due to both dark and stellar components. Remarkably, the maximum value of the circular velocity 
occurs at about 15 ± 3 kpc, independent of the galaxy mass: this seems to be a main kinematic imprint of the DM - luminous 
mass interaction occurring in spirals. Furthermore, Fig. 3 shows that the "Cosmic Conspiracy" paradigm has no observational 
support: there is no fine tuning between the dark and the stellar structural parameters to produce the same particular RC 
profile in all objects (e.g. a fiat one). Conversely, a number of relationships between the various structural parameters produce 
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Figure 4. The Universal Rotation Curve, normalized at its virial value Vurc {Rvir), as a function of normalized dark radius x = R/R v { r . 
Each curve, from the highest to the lowest, corresponds to M v i r defined as in Fig. 3. The bold line is the NFW velocity profile (see text). 

a variety of RC profiles. Moreover, the peak velocity of the stellar component Vj^* = Vd(2.2.Rd) = GMd/Rd k, with 
k = const, is not a constant fraction of the virial velocity as is found in ellipticals, (i.e a oc Vvir), but it ranges between the 
values 1 and 2 depending on the halo mass. 

Moreover, as in the NFW (and Burkert) RC profiles, the URC profiles are found (moderately) decreasing over most of 
the halo radial extent. The paradigm of flat rotation curves is obviously incorrect even/especially intended as an asymptotic 
behavior at large radii. In fact, we find that both U(0.05iiW),the velocity at the farthest radius with available kinematical in 
PSS and V(3Rd), a main reference velocity of the luminous regions of spirals, are significantly (10% — 30%) higher than the 
(observational) value of V V i r . This rules out a V = constant extrapolation of the inner RCs out to regions non mapped by 
the kinematics and DM dominated regions. We note that this result is independent of the adopted halo density profile and is 
far from being granted on theoretical grounds. 

In Fig. 4 we frame the URC from a full DM perspective by plotting Vjjrc{R/ Rvir', M„j r ). We set the virial mass M v i r 
as the galaxy identifier and R/R V i r as the radial "dark" coordinate, thus normalizing the amplitudes by V V i r oc M^. This 
ensemble of curves, a main goal of the present work, is parallel to those emerging in iV-body simulations and aims to represent 
the actual velocity profiles of spirals. In these variables the DM halos are self-similar; the whole system is self-similar in 
the outer regions, while in the innermost 30% of the halo size the baryons have influenced the dynamics and broken the 
self-similarity. In these coordinates it easily emerges that the maximum of the RC occurs at very different radii, viz. at ~ 2Rd 
for the most massive objects and at ~ 10-Rd for the least massive ones. Then, no reference circular velocities, to be considered 
as the actual physical counterparts of the empirical velocities of the Tully-Fisher relationship, exist in actual galaxy rotation 
curves. 

In Fig. 5 we zoom into the URC to look for the inner (luminous) regions of spirals from a baryonic perspective: the URC 
is so expressed as a function of the "baryonic" radial coordinate r/Ro- This figure corresponds to Fig. 4 of PSS, with the 
important difference that here the virial mass, rather than the galaxy luminosity, is the galaxy identifier. Plainly, an inverse 
correlation between the average steepness of the RC slope and the halo mass holds, similar to the slope-luminosity relationship 
found by Persic & Salucci (1988). In this coordinate the stellar matter is closely self-similar, and the different shapes of the 
RCs curves are mainly due to the M v i r — Md relation. In the space defined by normalized circular velocity - dark radius - 
halo mass, spirals do not occupy random positions, but a well defined plane of very small thickness. We clearly see that, by 
filling only less than 10 _3 % of the available volume, the available kinematics of spiral galaxies defines the Universal Rotation 
Curve. Let us notice that, in principle, theories of the formation of spirals do not trivially imply the existence of such a surface 
that underlies the occurrence of a strong dark - luminous coupling. 

We now show the URC DM density distribution. In Fig. 6 we show it as a function of x and M V i r . For x < 10 _1 
the well known core-cusp discrepancy emerges, i.e. the DM density of actual halos around spirals is about one order of 
magnitude smaller and radially much more constant than the NFW predictions. At x > 0.4, for a concentration parameter 
c = 13(M„ir/10 12 Mq)~ ' 13 , the observed halo densities are consistent with the NFW predictions for halos of the same virial 
mass. Note that this is a direct test: for halos with density profiles at x > 0.5 very different from the Burkert or the NFW 
profiles, Eq. 9b does not have solution. 

More specifically, let us constrain the analytical form of the outer DM distribution. For 2 < r/ro < 18, the following 
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Figure 6. The URC halo density vs. the NFW halo density of objects of the same mass, as a function of normalized radius and virial 
mass. The axes labels are x, logM v i r /Mq and log(p/(g cm -3 )) 
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approximation for the Burkert and NFW profile holds (y = r/ro, e = 0) 

2 06 ?; 0,86 

V URCH {y) = V URC H(3.24) i5Q+ y ylig+e . (11) 

Let us suppose that the actual outer DM velocity profile is different from the Burkert/NFW given by Eq. (11), i.e. e / 0. 
Then in Fig 7 we show that, even assuming large uncertainties in V v i r , in order to match both V(Ri) and V v i r , we must have 
e < 0.1. This is a first direct support for the Burkert and the NFW density law to be able to represent the outer regions 
(0.3 < r/R vir < 1) of DM galaxy halos. 

Notice that weak-lensing shear fields, at several hundreds kpc from the galaxy centers, are found compatible with the 
predictions of the NFW density profile, but cannot exclude non-NFW profiles ( Kleinheinrich ., et al., 2006 and references 
therein). 



4 DISCUSSION AND CONCLUSIONS 



In this paper, we have built the Universal Rotation Curve of spiral galaxies by means of kinematical and photometric data. 
We physically extended the URC, established for the inner region of galaxies in PSS out to R vir and have been able to employ 
the virial mass M vir as the parameter that characterizes spiral galaxies, and the virial radius R vir as a unit of measure for 
the radial coordinate. This URC is meant to be the observational counterpart of the NFW rotation curve, emerging from 
cosmological simulations performed in the CDM scenario. The URC yields the gravitational potential at any radius and it 
allows to link the local properties in the inner luminous regions with the global properties of the DM halos. 

DM halos have one (and likely just one) characteristic length scale, ro oc M®'£., which it is not naturally present in current 
scenarios of galaxy formation. Thus, they do not show any sign of an inner cuspy region of size r s oc Mj^ 4 . The halo velocity 
contribution Vurch rises with radius like a solid body at r ~ 0, decelerates to reach a maximum at 3.24 ro from where it 
start to slowly decrease out to R v i r with a slope that it is consistent with that of the NFW halos. The main significance of 
the URC concerns the full mass distribution (MD). First, it is possible to immediately exclude the following scenarios (and 
combinations of them): i) individual behavior, every object has its own MD; ii) unique behavior, every object has almost the 
same MD. Instead, the MD in spirals shows a remarkable mass-dependent systematics: both the dark and the stellar matter 
are distributed according to profiles that are functions of the total mass M vir (see Fig. 8). Finally, the DM halo becomes the 
dominant mass component in galaxies at different radii, according to the galaxy mass: from ~ W~ 2 R vir for the lowest masses, 
to ~ W~ 1 R V i r for the highest ones. 

We write the equilibrium velocity of the halos around spirals as the following approximation of the relations in the 
previous Sect.: 



Vurch = A(M vi . 



" 2 M 



1 + l(M vir 



tan 



1 



In 



1 + y(M v i r )x' 



with A(M vir ) = 0.406+1.08 Zo 5 [M„i r /(10 n Mo)]-0.688 {log[M vir /(10 11 M Q )]} 2 )+0.766 {log[M vir /(lO 11 M )]} :i and7(M vi 
26.78 [Mt, ir /(10 11 Mo)]~ ' 246 . This is the observational counterparts of JV-body outcomes. 

A Mathematica code for the figures in this paper is available at: http://www.novicosmo.org/salucci.asp. 



5 APPENDIX 

In this Appendix we discuss the observational evidence for the URC claim, the nature and the implications of which it is 
worth to clarify. The paradigm states that, when binned by luminosity, the RCs form a set of smooth, low-scatter synthetic 
curves, whose profiles and amplitudes are strong functions of the luminosity bin. 4 . Furthermore, the URC paradigm implies: 
i) rotation velocity slopes vs. rotation velocity amplitudes relationships (see fig (2) and (3) of PSS) and ii) a set of relations 
(Radial Tully Fisher relationship) holding at different radii x, defined as: 

log V(x) = a x M + b x 

where x = R/Rd and a x and b x are the fitting parameters, and M is the galaxy magnitude (Persic and Salucci 1991). 

Evidence for the URC claim and/or its above implications comes from: a) detailed analyzes of independent samples: 
Catinella et al 2006, (2200 RCs, see their Fig. 12), Swaters 1999, (60 extended RCs, see Chapter 4); b) independent analyzes 
of the PSS sample: Rhee (1996), Roscoe (1999); c) the finding of a very tight RTF in PSS and other three different samples 
Willick (1997, see below), Yegorova et al. (2007). 

The claim has been also tested by comparing the RCs of two samples of spirals (Courteau, 1998, 131 objects; Verheijen, 



4 The analytical form of the URC is built by assuming reasonable disk-halo velocity profiles, with three free parameters (V(Ropt), a, j3) 
that are obtained by x 2 -fitting the sinthctic curves. 
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Figure 7. Halo velocities at Ri and R v i r (filled circles) vs. the URC-halo velocity, given by Eqs. (11) (solid line) and vs. velocity profiles 
with average logarithmic slope steeper or shallower by an amount e = 0.1 (dashed line). 
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Figure 8. The dark halo (solid line), the disk (short dashed line), and the total (long dashed) URC mass profile for reference masses of 
IO^Mq and 10 13 Mq. The distribution for the intermediate halo masses can be derived from section 2 

1997, 30 objects) with the circular velocities predicted by the URCo, once that the values of galaxy luminosity and disk length- 
scale are inserted in it. The face- value result of the test: 2/3 of the RCs are in pretty good agreement with the universal 
curve, while 1/3 show some disagreement, indicates that the URC is a useful tool to investigate the systematics of the mass 
distribution in spirals, but also it questions about its universality. However, while some of this disagreement may reflect an 
inefficiency of the URCo to reproduce the RCs, the actual performance of the URC is better than it is claimed. In fact, 
spurious data vs predictions disagreements are created in performing this test and precisely when they insert in URCo the 



